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Abstract 

We present exact calculations of the Potts model partition function Z{G, q, v) 
for arbitrary q and temperature-like variable v on n- vertex strip graphs G of the 
honeycomb lattice for a variety of transverse widths equal to Ly vertices and 
for arbitrarily great length, with free longitudinal boundary conditions and free 
and periodic transverse boundary conditions. These partition functions have 
the form Z{G,q,v) = Ylj=i^'^ ^z,G,j{^z,G,j)^ ■, where m denotes the number of 
repeated subgraphs in the longitudinal direction. We give general formulas for 
Nz,G,j for arbitrary Ly. We also present plots of zeros of the partition function 
in the q plane for various values of v and in the v plane for various values of q. 
Explicit results for partition functions are given in the text for Ly = 2,3 (free) 
and Ly = A (cylindrical), and plots of partition function zeros are given for Ly 
up to 5 (free) and Ly = 6 (cylindrical) . Plots of the internal energy and specific 
heat per site for infinite-length strips are also presented. 
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1 Introduction 



In this paper we present some theorems on structural properties of Potts model 
partition functions on strips of the honeycomb (he) lattice of arbitrary width equal 
to Ly vertices and arbitrarily great length. We also report exact calculations of these 
partition functions for a number of honeycomb-lattice strips of various widths and 
arbitrarily great lengths. Using these results, we consider the limit of infinite length. 
For this limit we calculate thermodynamic functions and determine the loci in the 
complex q and temperature planes where the free energy is non-analytic. This is an 
extension of our previous study for this lattice in [1]. 

We briefly review the definition of the model and relevant notation. Consider 
a graph G = {V,E), defined by its vertex set V and edge set E. (Here we keep 
the discussion general; shortly, we will specialize to strips of the honeycomb lattice.) 
Denote the number of vertices and edges as = n and \E\, respectively. On the 
graph G, at temperature T, the Potts model is defined by the partition function 
Z{G,q,v) = ^{o-^} e"^"^ with the (zero-field) Hamiltonian H = —JJ2{ij) where 
cTj = 1, . . . ,g are the spin variables on each vertex i E V; P = (ksT)^^; (ij) G E 
denotes pairs of adjacent vertices, and J is the spin-spin interaction constant. We 
use the notation K = pj, a = , and v = — 1. The physical ranges are thus (i) 
a > 1, i.e., f > corresponding to oo > T > for the Potts ferromagnet with J > 0, 
and (ii) < a < 1, i.e., —l<v<0, corresponding to < T < oo for the Potts 
antiferromagnet with J < 0. Let G' = {V,E') with E' C E. Then Z{G,q,v) can be 
defined for arbitrary q and v by the formula p] 

z{G,q,v)= J2 g'^'^'^t^i'^'i (1.1) 

G'CG 

where k{G') denotes the number of connected components of G'. We define a (re- 
duced) free energy per site / = —^F, where F is the free energy, via 

f{{G},q,v)= limln[Z(G,g,t;)V"] , (1.2) 

n^oo 

where the symbol {G} denotes lim^^oo G for a given family of graphs G. 

Our exact results on Z{G,q,v) apply for arbitrary q and v. We consider free 
and cylindrical strip graphs G of the honeycomb lattice of width Ly vertices and 
of arbitrarily great length vertices. Here, free boundary conditions (sometimes 
denoted FF), mean free in both the transverse and longitudinal directions (the latter 
being the one that is varied for a fixed width), while cylindrical boundary conditions 
(sometimes denoted PF) mean periodic in the transverse direction and free in the 
longitudinal direction. We represent the strip of the honeycomb lattice in the form 
of bricks oriented horizontally. For the honeycomb lattice with cylindrical boundary 
conditions, the number of vertices in the transverse direction, Ly, must be an even 
number, and the smallest value without degeneracy (multiple edges) is Ly = 4. Exact 
partition functions for arbitrary q and v have previously been presented for strips of 
the honeycomb lattice with free boundary conditions for width Ly = 2 in Our 
new results include theorems that describe the structure of the Potts model partition 
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function for strips with free and cylindrical boundary conditions, of arbitrary width 
and length and explicit calculations using the transfer matrix method (in the Fortuin- 
Kasteleyn representation f3^) for strips with free and cylindrical boundary conditions 
with width Ly = 3 (free) and Ly = 4 (cylindrical). We have carried out similar 
calculations for Ly < 7 (free) and Ly = 6 (cylindrical); these are too lengthy to 
include here. We shall also present plots of partition function zeros in the limit of 
infinite length, for widths 2 < Ly < 5 (free) and Ly = 4,6 (cylindrical). Related 
calculations of Potts model partition functions for arbitrary q and v on fixed-width, 
arbitrary-length strips of the square and triangular lattices are [l]-[ini and [TU t fTT I fT^. 
respectively. Analogous partition function calculations for arbitrary q and v on finite 
sections of 2D lattices with fixed width and length include flbi [T6] . The special 
case f = — 1 is the zero-temperature limit of the Potts antiferromagnet, for which 
Z{G,q, —1) = P{G,q), where P{G,q) is the chromatic polynomial expressing the 
number of ways of coloring the vertices of the graph G with q colors such that no two 
adjacent vertices have the same color. 

As part of our work, we calculate zeros of the partition function in the q plane for 
fixed V and in the v plane for fixed q. In the limit of infinite strip length, L^ — »• oo, 
there is a merging of such zeros to form continuous loci of points where the free 
energy is nonanalytic, which we denote generically as B. For the limit L^ — ^> oo of a 
given family of strip graphs, this locus is determined as the solution to an algebraic 
equation and is hence an algebraic curve. 

There are several motivations for this work. Clearly, exact calculations of Potts 
model partition functions with arbitrary q and v are of value in their own right. This is 
especially true since there are no exact calculations of the free energy f{{G}, q, v) for 
arbitrary q and v on an infinite lattice of dimension two or higher. Exact calculations 
on lattice strips of fixed width and arbitrarily great length thus provide a useful set 
of results complementing other methods of analysis such as series expansions and 
Monte Carlo simulations of the Potts model. Our structural theorems elucidate the 
form of the partition function on these strips for arbitrarily great widths as well as 
lengths. The honeycomb lattice is of interest since, together with the square and 
triangular lattices, it comprises the third and last regular tiling of the plane which is 
homopolygonal, i.e. composed of of a single type of regular polygon. While critical 
properties describing the second-order phase transition of the Potts ferromagnet are 
universal and independent of lattice type, the behavior of the Potts antiferromagnet 
is sensitively dependent on lattice type, so that studies of this model on different 
lattices and lattice strips are valuable. 

There is a particular motivation for carrying out exact calculations of Potts model 
partition functions with arbitrary q and f , because this allows one to investigate more 
deeply a unique feature of the model, which is qualitatively different from the behav- 
ior on either the square or triangular lattice, namely the property that the critical 
temperature of the Potts antiferromagnet on the honeycomb lattice decreases to zero, 
i.e., the critical v decreases to —1, at a non-integral value, q = (3 + v^)/2 (as reviewed 
below in connection with the criticality condition, eq. (16.31) ). This is formal, since 
for q ^ Z_|_, the model (with either sign of J) is only defined via the representation 
(11. ip . and in the antiferromagnet case (i.e., for — 1 < f < 0), this formula can yield 
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a negative, and hence unphysical, result for the partition function. One obviously 
cannot investigate this formal criticality using the Hamiltonian formulation, which 
requires g G Z+. 

Moreover, calculations of complex-temperature zeros of the partition function 
show how the physical phases can be generalized to regions in the plane of a complex- 
temperature variable, as was discussed for the 2D Ising model on the square lattice 
p!3| [Ti] . Calculations of partition function zeros on long finite lattice strips, and the 
loci B in the infinite-length limit, also yield interesting insights into properties of the 
corresponding phase diagrams in the complex-temperature and complex-g planes. 



2 General Structural Theorems 

2.1 Preliminaries 

In this section we prove several general theorems that describe the structure of 
the partition function for the honeycomb-lattice strips under consideration. Let m 
denote the number of bricks in the longitudinal direction for such a strip. Then the 
length (number of vertices in the longitudinal direction) is 



2m + 1 for odd L^. 



(2.1) 



2m + 2 for even L^. . 
For this type of strip graph, Z{G, q, v) has the form 

Z{G,q,v)= ^G,{\z,G,r 
i=i 



(2.2) 



where the coefficients cgj and corresponding terms Xgji as well as the total number 
Nz,G,\ of these terms, depend on the type of strip (width and boundary conditions) 
but not on its length. In the special case t> = — 1, the numbers Nz,g,\ will be denoted 
Np,G,\- We define Nz,hc,BCy BCx,Ly,\ as the total number of A's for the honeycomb- 
lattice strip with the transverse and longitudinal boundary conditions BCy and BC^ 
of width Ly. Henceforth where no confusion will result, we shall suppress the A 
subscript. The explicit labels are Nz,hc,FF,Ly and Nz,hc,PF,Ly for the strips of the 
honeycomb lattices with free and cylindrical boundary conditions. 



2.2 Case of Free Boundary Conditions 

Theorem 2.1 For arbitrary Ly, 

Ciy for odd Ly , 



N. 



ZM,FF,L 



y 



^Ly + (4/2) /^'^ e^s'^ 



(2.3) 
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Proof A honeycomb-lattice strip with free boundary conditions is symmetric 
under reflection about the longitudinal axis if and only if Ly is even. Therefore, for 
odd Ly, the total number of A's in the Potts model partition function, Nz,hc,FF,Ly, 
is the same as the number for the triangular lattice strips Nz^tri,FF,Ly, namely, the 
number of non-crossing partitions of the set {1, 2, Ly}. This is the Catalan number 
[T9l[20], Cl = (I/-|-l)~^(^^). For even L^, the reflection symmetry reduces Nz,hc,FF,Ly 
to the number for the square lattice strips Nz,sq,FF,Ly, which was given in Theorem 5 
of [7]. □ We list the first few values of Nz,hc,FF,Ly in Table [H 



Table 1: Numbers of A's for the Potts model partition function and chromatic polynomials for the 
strips of the honeycomb lattices having free boundary conditions and various widths Ly. 



Ly 


Nz,hc,PF,Ly 


Np,hc,FF,Ly 


2 


2 


1 


3 


5 


3 


4 


10 


5 


5 


42 


19 


6 


76 


25 


7 


429 


145 


8 


750 


194 


9 


4862 


1230 


10 


8524 


1590 


11 


58786 


11139 


12 


104468 


14681 



We next discuss some combinatorics which will be used in our next theorem. 
Let us denote the elements in the fc'th column of the j'th row of Pascal's triangle 
as P{j,k), with the value P{j,k) = (^) (the binomial coefficient). The relation 
between these elements is P{j, k) = P{j — 1, k — 1) + P{j — 1, k) with P(0, 0) = 1, 
P{j, —1) = 0, and P{j, k) = for j < k. That is, each element is the sum of the two 
numbers immediately above it. For the reader's convenience, we display the first few 
rows of Pascal's triangle in Fig. [TJ Now, P{j, k) is the number of ways to place k 
black beads and j — k white beads in a line. Similarly, what is known as Losanitsch's 
triangle [T7] is given by the number of ways to put j beads of two colors in a line, 
modulo refiection symmetry. We denote the entry in the fc'th column of the j'th row 
of this triangle as L{j, k). For reference, the sequence formed by reading the entries 
in this triangle (from left to right) by rows is listed as sequence A034851 in [T7]. The 
relation between these entries is essentially the same as for Pascal's triangle, except 
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when j is even and k is odd: 

k) = L{j -l,k-l)+ L{j -l,k)- 6j even.fc odd (^^^''_^ ^^^^2) (2-4) 

where Sj even.fe odd = 1 if j is even and k is odd and zero otherwise. For reference, the 
first few rows of Losanitsch's triangle are shown in Fig. [2l Now form a new triangle 
by subtracting the entries in Losanitsch's triangle from the corresponding entries in 
Pascal's triangle, and denote its elements as PL{j,k). The first few rows of this 
triangle are displayed in Fig. [31 
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1 8 28 56 70 56 28 8 1 

Figure 1: Pascal's triangle. 
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14 12 19 19 12 4 1 

1 4 16 28 38 28 16 4 1 

Figure 2: Losanitsch's triangle. 



The next theorem concerns the number of A's Np^hc,FF,Ly in the chromatic poly- 
nomial for the free he strip. 
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Figure 3: Triangle formed by subtraction of elements of Losanitsch's triangle from corresponding 
elements of Pascal's triangle. 



Theorem 2.2 For arbitrary Ly, 



y ' 



E£r'"'A/,._,_.(<'.-'>/^) for odd L 

Ei^=f ^ PL{Ly/2 - 1, i)Np^FP,i(Ly+i-i)/2] for even Ly . 



where Ml is the Motzkin number IT^ . Np^sq,FF,Ly is the number of X 's for the square 
lattice with free boundary conditions given in Theorem 2 of fT], Np^FP,Ly is the total 
number of A 's for the square, triangular, or honeycomb lattices with cyclic boundary 
conditions given in eq.(5.2) of \2D^ . and PL{j,k) is the number given by the sub- 
traction of the elements of Losanitsch 's triangle from the corresponding elements of 
Pascal's triangle. 

Proof Consider odd strips widths Ly first; for these, there is no refiection sym- 
metry. For each transverse shce containing Ly vertices of the honeycomb lattice with 
free boundary conditions, there are {Ly — l)/2 edges. Compared with the calculation 
of partitions for strips of the square lattice, the same numbers of edges are removed, 
such that the end vertices of each of these missing edges are allowed to have the same 
color. Firstly, all of the possible partitions for the triangular-lattice strip with free 
boundary conditions Np^tri,FF,Ly = ^Ly-i [IS] are valid for the honeycomb lattice. 
Among {Ly — l)/2 missing edges, if one pair of end vertices has the same color, the 
number of partitions is given by Np^tri,FF,Ly-i- As an example, for the Ly = 3 slice 
of the honeycomb-lattice strip, there is an edge connecting vertices 1 and 2 and no 
edge connecting vertices 2 and 3. It follows that the set of partitions is comprised of 
{1,61^3,62^3} in the shorthand notation used in [71 |T2]. Similarly, if two pairs of end 
vertices of missing edges separately have the same color, there are [(•^y-^y^) choices 
for the locations of missing edges and the number of partitions for each choice is 
Np tri,FF,Ly-2- By including all the possible missing edges with the same color on 
their end vertices, the first line in eq. (12.51) is established for the honeycomb lattice 
with odd Ly. 
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For the honeycomb lattice with even Ly, reflection symmetry must be taken into 
account. We thus consider partitions with the transverse slice composed of edges 
connecting vertices 1 and 2, vertices 3 and 4, vertices Ly — 1 and Ly. That is, 
there are Ly/2 edges in each slice and Ly/2 — 1 missing edges. Naively, one would 

expect that the number of partitions is J2i=o^ ^ Np^sq,FF,Ly-i{^^^f~^) by the above 
argument. However, this would actually involve an overcounting, because certain 
partitions with the same color on end vertices of missing edge(s) are equivalent under 
the reflection symmetry. For example, the Ly = 6 strip of the honeycomb lattice 
has two missing edges, i.e., there is no edge between vertices 2 and 3, and between 
vertices 4 and 5. There are Np^sq,FF,5 = 7 partitions if vertices 2 and 3 have the same 
color, and separately seven partitions if vertices 4 and 5 have the same color. Five 
partitions for the first case ^2,3, ^2,3^1,65 ^2,3,5, ^2,3^1,4,6 and 62,3,561^ are equivalent 
under reflection to the following partitions for the second case: 64^5, S^^^Si^q, ^2,4,55 
^4,5^1,3,6, and 52,4,5^1,6- In general, the double-counted partitions are those symmetric 
partitions in Np^sq,FF,Ly-i such that i pairs of end vertices of missing edges have 
the same color. This number of symmetric partitions is |-/V'p,FP,[(La+i-i)/2]; given as 
Theorem 1 in [7] . The double counting occurs when two choices of i missing edges are 
reflection-symmetric to each other. This number of asymmetric choices of missing 
edges is given by PL{j,k). □ We hst the first few values of Np^hc,FF,Ly for the 
honeycomb lattice with free boundary conditions in Table [H 

2.3 Case of Cylindrical Boundary Conditions 

For the honeycomb lattice with cylindrical boundary condition, only strips with 
even Ly can be defined. The number of A's can be reduced from the number of non- 
crossing partitions Nz,tri,FF,Ly = ^Ly with length- two rotational symmetry, and this 
first reduction will be denoted as N'^^i^c,PF,Ly- This number can be further reduced with 
reflection symmetry, and this further- reduced number will be denoted as Nz,hc,PF,Ly- 
We list the first few relevant numbers in the Potts model partition function for the 
strips of the honeycomb lattices with cylindrical boundary conditions in Table El 



Lemma 2.1 For arbitrary even Ly, 

^A'Wf,.. - = J: *(^)('') (2-6) 



even d\Ly; l<d<L 



y 



where d\Ly means that d divides Ly, and (t){n) is the Euler totient function, equal 
to the number of positive integers not exceeding the positive integer n and relatively 
prime to n. 

Proof Because the honeycomb-lattice strip has length-two rotational symmetry, 
the value ^N'z^ 

hcPFLy Klust be larger than Nz,tri,FF,Ly- For strips with cylindrical 
boundary conditions, partitions can be classified according to the periodicity d modulo 
rotations. That is, a partition could be transformed back to itself when the Ly vertices 
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Table 2: Numbers of A's in the Potts model partition function for the strips of the honeycomb 
lattices having cylindrical boundary conditions and various even Ly. 



Ly 


Nz,tri,FF,Ly 


^^Z,hc,PF,Ly 


Nz,hc,PF,Ly 


'^Nz,hc,PF,Ly 


2 ^^ZM,PF,Ly 
—Nz,tri,FF,Ly 


2 


2 


2 


2 


2 





4 


14 


10 


8 


6 


6 


6 


132 


48 


34 


20 


12 


8 


1430 


378 


224 


70 


82 


10 


16796 


3364 


1808 


252 


24 


12 


208012 


34848 


17886 


924 


1076 



are rotated by length where d denotes any of the positive integers that divide Ly. 
The number of partitions which have periodicity d modulo these rotations was denoted 
as 2ad in the proof of Theorem 2.2 of [T5], and is given by 




where is the Mobius function, defined as ^{n) = — 1 if is prime, /i(n) = if n 
has a square factor, and /i(n) = 1 for other n. Now the periodicity d can be either 
odd or even. For a partition with odd periodicity rf, the contribution to the excess of 
{Ly/2)N'zi^^ppj^^ relative to Nz^tri,FF,Ly is 2ad{Ly/2 — d). On the honeycomb-lattice 
strip with length-two rotational symmetry, a partition with even periodicity d and its 
length-one rotation are not equivalent. Thus for a partition with even periodicity d 
and its length-one rotated counterpart, the contribution is 2ad{Ly — d). We have 

Y^'zM,PF,Ly-Nz,tri,FF,Ly = '^MLy - d) + ^ 2ad{^ - d) 

even d\Ly odd d\Ly 

= ^2arf(Ly-(i)- ^ UdLy 

d\Ly odd d\Ly 

J2 HLy/d)(^f\- ^'^^y (2.8) 

d\Ly; l<d<Ly ^ ^ odd d\Ly 

where the last line follows from the proof of Theorem 2.2 in Using eq. fl2.7p for 
the second summation in eq. fl2.8p is 

odd d\Ly 



odd d\Ly d'\d ^ ' 

E E "^1^)0 



odd d'\Ly odd d\Ly;d'\d 
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= E lO E §M</"). (2.9) 

odd d'\Ly ^ ^ odd d"\L'y 

where we change the variables to c?" = djd^ and L'^ = Ly/d'. (This is given as sequence 
A062570 in [T7].) The summation of d" is 

J2 = 0(2^^;) = 2<p{Ly/d') , (2.10) 

odd d"\L'y 

where that second identity follows because L'y is even by the formula for the Euler 
totient function, (f){n) = ''^^primo p\n(^ ~ Vp)- Therefore, we find 

c^dLy= Yl 'PiLy/d)!^^), (2.11) 

odd d\Ly odd d\Ly 

and the proof is completed. □ 



Lemma 2.2 For arbitrary even Ly, 

'2Nz,hc,PF,Ly - Nz^f^c,PF,Ly = Nz,FP,^ (2-12) 

where Nz^FP,Ly = {^i^) the total number of X's for the square or triangular or 
honeycomb lattices with cyclic boundary conditions, given in eq. (5.6) of f2Dj . 

Proof We use the result in Theorem 2.1 of [12] that 2Nz,sq,PF,Ly — Nz,tri,PF,Ly = 
N^ppLy_ for even Ly is the number of partitions with both rotation and reflection 

symmetries. The partitions that do not have reflection symmetry appear once in 
Nz,sq,PF,Ly but twice (itself plus its mirror image) in Nz,tri,PF,Ly, so that they do not 
contribute to 2Nz,sq,PF,Ly — Nz,tri,PF,Ly- It was shown there that for even L^, there 
are two classes of reflection symmetries, denoted as type I and type II. For type I 
partitions, the reflection axis does not go through any vertex; while for type II parti- 
tions, the reflection axis goes through two vertices. There are at least two partitions 
that belong to both type I and II classes, namely, the partitions 1 (identity) and 
^1,2,. ..,Ly (unique block). Denote the set of partitions of Nz^tri,PF,Ly as Pz,tri,PF,Ly, 
that of N'zhcPFL,, as P'zhcPFL,,^ and that of Nz,hc,PF,Ly as Pz,/ic,pf,l„- For exam- 
pie, Pz,tri,PFA = {l,5i,2,3,45'^i,2,5i,3, 51,2^3,45 '^i,3,4}- As the Cylindrical strip of the 
honeycomb lattice only has length-two rotation symmetry, N'^ pp is larger than 
Nz,tri,PF,Ly- The excess partitions in P'z,hc,PF,Ly relative to P z,tri,PF,Ly can be obtained 
by making length-one rotations for certain partitions in P z,tri,PF,L ■ For example, in 
addition to those in P z,tri,PF,Ly^ '^'zhcPFi contains ^2,35 ^2,4, '^2,35i,45 ^1,2,4- Notice that 
the partitions that belong to both type I and II classes remain the same after the 
length-one rotation (modulo length-two rotational symmetry), so that they should be 
excluded in the doubling. For the honeycomb lattice, the reduction from P'zhcPfl 
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to Pz,hc,PF,Ly occurs Only for type I partitions. Considering the Ly = 4 strip again 
as an example, we observe that 61^3 and ^2,4 are equivalent, so that only one of them 
should be kept for PzM,PF',Ly, and similarly for ^i^s^ and 5i,2,4- Since the numbers of 
partitions in type I and II classes are the same, the right hand side of eq. fl2.12l) is 
the same as that in Theorem 2.1 of |12| for even Ly. □ 



The exact formula for Nz,hc,PF,Ly follows from Lemmas 12. II and 12. 2t 
Theorem 2.3 For arbitrary even Ly, 



Z,hc,PF,Ly 



2 \Ly/2j Ly 



even d\Ly; l<d<Ly 



Ly 

d 



2d 

d 



(2.13) 



We list the first few relevant numbers in the chromatic polynomial for the strips of 
the honeycomb lattices with cylindrical boundary conditions in Table [31 Analogously 
to Lemmas 12. II and I^T^ we state the following conjectures. 



Table 3: Numbers of A's for the chromatic polynomial for the strips of the honeycomb lattices 
having cylindrical boundary conditions and various even Ly. 



Ly 


np{hc, Ly, 0) 


^^P,hc,PF,Ly 


Np^hc,PF,Ly 


'^Np^hc,PF,Ly 
^^P,hc,PF,Ly 


hi AT' 

2 ^^P,hc,PF,Ly 

—np{hc, Ly, 0) 


2 


1 


1 


1 


1 





4 


6 


5 


4 


3 


4 


6 


43 


17 


12 


7 


8 


8 


352 


99 


62 


25 


44 


10 


3114 


626 


346 


66 


16 


12 


29004 


4907 


2576 


245 


438 



Conjecture 2.1 For arbitrary even Ly, 
Y^P,hc,PF,Ly - np{hc, Ly, 0) = 



^ ) Np^r^cm (2-14) 



where np{hc, Ly, 0) is the number of A 's for the cyclic strips of the honeycomb lattice 
with level d = 0, and NphcFPL is the total number of \ 's for these strips given in 



Compared with Lemma [2711 the number of non-crossing partitions Nz^tri,FF,Ly = Ciy 
is now replaced by the corresponding number np{hc, Ly,0) for the honeycomb strip in 
the chromatic polynomial, and the number (^^) = Nz^hc,FP,d is replaced by Np^hc,FP,d- 
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Conjecture 2.2 For arbitrary even Ly, 

\^Phr FP^ /'^^ Ly/2> 1 , 



'2Np^hc,PF,Ly - Nphcppr 



for even Ly/2 . 

(2.15) 



Conjectures 12.11 and 12.21 imply an exact formula for Np^hc,PF,Ly- 

3 Potts Model Partition Functions for Strips of 
the Honeycomb Lattice with Free Boundary Con- 
ditions 

The Potts model partition function for a strip of the honeycomb lattice of width 
Ly and length vertices with free boundary conditions is given by 

Z{Ly X L,, FF, g, v) = w^-J"'- u,^ (3.1) 

where T = V ■ H2 ■ V ■ Hi is the transfer matrix. Hi and H2 are matrices corresponding 
respectively to adding two kinds of transverse bonds in a slice, and V corresponding 
to adding longitudinal bonds in each slice. The number m is related to as defined 
in eq. (12.11) . and the vector w is given by 

{wJdd = v"^ ■ Hi for odd L^, , 

(3.2) 
wJven = v'^ ■ H2 ■ V ■ Hi for evcu . 

Hereafter we shall follow the notation and the computational methods of [71 [T2] . 

The matrices T, V, Hi and H2 act on the space of connectivities of sites on the 
first slice, whose basis elements v-p are indexed by partitions V of the vertex set 
{1, . . . , Ly}. In particular, Ujd = ^{{i},{2},...,{Ly}}- We denote the set of basis elements 
for a given strip as P = {v-p}. 

An equivalent way to present a general formula for the partition function is via a 
generating function. Labelling a lattice strip of a given type and width as Gm, with 
m the length, one has 

00 

r(G, q, v,z) = Y, z^Z{G^, q, v) (3.3) 



m=0 

where T{G, q, v, z) is a rational function 

r(G, g, f , z) = — — (3.4) 

V{G,q,v,z) 

with 

degJAT) 

AfiG,q,v,z)= J2 ^G,jZ^ (3.5) 

j=0 
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Nz.hc,BC.Ly Nz,hc,BC.Ly 

V{G,q,v,z) = l+ Yl bG,jz'= n i'^-^z,G,jz) (3.6) 

j=i j=i 

where the subscript BC denotes the boundary conditions. In the transfer-matrix 
formahsm, the A^.gj's in the denominator of the generating function, eq. (13. 6p . are 
the eigenvalues of T. 

Strips of the honeycomb lattice with free boundary conditions are well-defined 
for widths Ly > 2. The partition function Z{G,q,v) was calculated, for arbitrary 
q, V, and m, for the strip with Ly = 2 and free boundary conditions in [T], using 
a systematic iterative application of the deletion-contraction theorem. Here after 
re-expressing the results for Ly = 2 in the present transfer matrix formalism, we 
shall report explicit results for the partition function for strips with Ly = 3 and free 
boundary conditions. For 4: < Ly < 7, the expressions for T{Ly), w{Ly) and Ui^^^Ly) 
are too lengthy to include here. They are available from the authors on request. 



3.1 Ly — 2 

For the strip with width Ly = 2, we only have to consider odd L^. The number 
of elements in the basis is equal to C2 = 2: P = {l,5i^2}- In this basis, the transfer 
matrix and the other relevant quantities are given by 



T 



w 



odd 



Rll D1F2R12 



(3.7a) 
(3.7b) 



where 



Dk = V + k 
Fk = q + kv 

R12 = q^ + 2qv + 2v^ . 



(3.8a) 
(3.8b) 
(3.8c) 
(3.8d) 



In terms of this transfer matrix and these vectors, one calculates the partition function 
Z{Gm,q,v) for the strip with a given length m via eq. (13.11) . Equivalently, one can 
calculate the partition function using a generating function, and this was the way in 
which the results were presented in [Ij, with V = n?=i(l ~ ^hcf2,jz) and 



where 



and 



A 



/ic/2,(l,2) 



Ml ± a/Ms 



Ml 



+ 5q^v + 10q\^ + Wqv^ + + 



(3.9) 
(3.10) 



Mo 



q^ + lOq'^v + 45g%2 + 120q^v^ + 208gS^ - 2g%^ + 2AAqh^ - Qqh^ 
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+196gV - 4g2t;7 + IQAqv^ + Aqv^ + 32v^ + 8v^ + . (3.11) 

The product of these eigenvalues, i.e., the determinant of T, is 

det(T) = v\l + v){v + qf = v^D^F^ . (3.12) 

The vanishing of this determinant at f = — 1 and v = —q occurs because in each 
case one of the two eigenvalues is absent for, respectively, the chromatic and flow 
polynomials [211 122]- Analogous formulas can be given for det(T) for higher values of 
Ly] we omit these for brevity. 



3.2 



Ly — 3 



For the he strip of width Ly = 3 the number of elements in the basis 
merating partitions is given by the Catalan number C3 = 5. This basis 
{1, 5i^25 (^1,35 (^2,35 '5i,2,3}- In this basis, the transfer matrices and the other 
quantities are 



for enu- 
is P = 
relevant 



T 



FiF2S'ii 


DiF.Su 




Su 




v^FiF2 


v^DiFiF2 


V'S23 


V'S2, 


v^D,S2, 


v^Fi 


v^DiFi 




V'S2, 


v^D,S2, 


v^FiF2Ri2 


v^DiFiS42 


V^F2S,3 






v'F^ 


v^DiFi 


V'S2, 







\ 



T 

Wodd 



T 
Uid 



v 

q {q{q + v),q{l + v),q + v,q + v,l + v) 
q [Fl D,F^, FiXi, F^Xa, 
(1,0,0,0,0) 



(3.13a) 

(3.13b) 
(3.13c) 
(3.13d) 



where the Sa and Xi. are defined in a shorthand notation as 



5*11 


= + 5q^v + llq^v^ + 12qv^ + 7v'^ 


(3.14a) 


>S'l2 


= + 6q\ + 15q^v^ + 19gw^ + llw^ 


(3.14b) 


>S'l3 


= + 4gt> + 5f ^ 


(3.14c) 


o> 
^13 


= q^ + Aq^v + Iqv"^ + + 


(3.14d) 


5*14 


= q^ + 8q\ + 28q\^ + q^v^ + bAq'^v^ + hq^v^ + 59gw^ + Qqv^ 


+ 32i;^ + 7v^ 






(3.14e) 


5'l5 


= q^ + 7q\ + 21q\^ + q\'^ + 33qv^ + Aqv^ + 2Av^ + 


(3.14f) 


>S'23 


= 2q + ^v + v'^ 


(3.14g) 


'S'25 


= g + 4t> + 


(3.14h) 


'S'33 


= g^ + 4gt; + Gv"^ + 


(3.14i) 


>S'42 


= g^ + 3gf + 3f ^ 


(3.14j) 


'S'43 


= g^ + 3gt; + + 


(3.14k) 


5*44 


= q^ + 6qh + q%^ + 12qv^ + 3qv^ + lOv^ + 3w^ 


(3.141) 


'S'45 


= g^ + 5qv + qv^ + 7f ^ + 2v^ 


(3.14m) 
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Xi = + Aq^v + Qqv"^ + Av^ + 
X2 = q^ + 3qv + 3v^ + 



(3.15a) 
(3.15b) 



We will discuss properties of the resultant partition functions below. 

4 Potts Model Partition Functions for Honeycomb- 
lattice Strips with Cylindrical Boundary Condi- 
tions 

For the honeycomb lattice with cylindrical boundary conditions, the width must 
be even (and larger than two in order to avoid the degenerate situation of vertical 
edges forming emanating from and returning to a given vertex). The Potts model 
partition function Z{G,q,v) for a honeycomb-lattice strip with cylindrical boundary 
conditions can be written in the same form as in eg. (13. II) . Here either Hi or H2 should 
include the bond connecting the boundary sites in the transverse direction. The 
dimension of the transfer matrix can be reduced by the two symmetries discussed in 
Section [21 namely the length- two translation symmetry along the transverse direction 
and reflection symmetry. This number is Nz,hc,PF,Ly and is given in terms of Ly 
by eq. fl2.13p (see Table [2] for some numerical values). We consider the basis in 
the translation-invariant and reflection-invariant subspace to construct the transfer 
matrix and the corresponding vectors. To simplify the notation, we will still use T, 
w and as in eq. fl3.ip . so that the partition function is given by the analogous 
equation, Z{Ly x L^;, P-F, g, v) = w"^ ■ T™ ■ Ui^. We have calculated the transfer matrix 
T{Ly) and the vectors w{Ly) and Uid{Ly) for Ly = A and Ly = 6. The explicit results 
for Ly = A are given in the appendix; the results for Ly = 6 are too lengthy to present 
here, and are available upon request. 

5 Partition Function Zeros in the q Plane 

In this section we shall present results for zeros in the g-plane for the partition 
function of the Potts antiferromagnet on strips of the honeycomb lattice with free and 
cylindrical boundary conditions, for various values of the temperature-like variable v. 
Fig. m shows these zeros for strips of widths 2 < L < 5 and free boundary conditions. 
In the limit Lx —>■ 00 the zeros merge to form sets of curves which, together, comprise 
the locus B. As an illustration of this, in Fig. [5] we show these zeros for the free 
strips of width Ly = 3, together with the loci B for various values of v. One sees that 
the strip lengths that we use to calculate the zeros are sufficiently great that most 
of these zeros lie rather close to the infinite-length asymptotic loci. This behavior 
- that partition function zeros calculated for long strips generally lie close to the 
asymptotic loci B - is similar to what we found in our previous work. Hence, one can 
draw a reasonably good inference for many of the features of these loci B from the 
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zeros. The corresponding partition-function zeros for honeycomb-lattice strips with 
Ly = 4,6 and cyhndrical boundary conditions are shown in Figure [61 Again, one 
could calculate the asymptotic loci B, but since the zeros already give a reasonably 
good idea of the structure of these loci, they will suffice for our present purposes. 
Our Figs, m and [6] include calculations up to Ly=5 and 6, respectively. We previously 
studied the case Lj^ = 2 in fl] (again for arbitrary q and temperature). From our 
present results, we see that, for general values of v, as the width Ly increases, the 
curve envelope moves outward somewhat and the arc endpoints on the left move 
slowly toward q = 0. This behavior is consistent with the hypotheses that for a given 
V, as Ly oo, (i) Bq would approach a limiting locus as Ly — > oo and (ii) this locus 
would separate the q plane into different regions, with a curve passing through g = 
as well as a maximal real value, qdv). This is qualitatively similar to the behavior 
that was found earlier for the square-lattice strips [U El E] and the triangular-lattice 
strips [12]. As expected, the convergence to the limit Ly ^ oo seems to be faster 
with cylindrical boundary conditions, as there are no surface effects when the length 
is made infinite. In this limit Ly —>■ oo, one expects that the locus B will have the 
property that the maximal point at which it crosses the real axis, qdv), is equal to the 
solution, eq. (16.91) (with the plus sign for the square root), of the criticality condition 
on the infinite honeycomb lattice, eq. (16. 3p below. 

We recall some results for the special case v = —1 corresponding to the zero- 
temperature Potts antiferromagnet. Chromatic zeros were calculated for a finite patch 
of the honeycomb lattice with free and cylindrical boundary conditions in (Fig. 8 of) 
[23] . Chromatic zeros and asymptotic loci for strips of this lattice were calculated 
for free longitudinal boundary conditions in [H [211 f*^^ periodic longitudinal 

boundary conditions in [H [201 [251 [2S] (see, e.g.. Fig. 6 of [21] for Ly = 3 and free b.c. 
Fig. 1 of [25] for Ly = 2 and cychc b.c. Fig. 17 of [I] and Fig. 7 of [26] for Ly = 3 
and cyclic b.c, and Figs. 8 and 9 of [26] for Ly = 4 and Ly = 5 with cychc b.c). 
Comparisons of the arcs on B for v = —1 obtained for strips with free longitudinal 
boundary conditions led to the inference that in the limit Ly ^ oo these loci would 
separate the q plane into regions including a curve passing through q = [211 ■ The 
property that, for v = —1, Bq separates the q plane into regions with one of the curves 
on Bq passing through the origin is also observed for lattice strips with finite width Ly 
if one imposes periodic longitudinal boundary conditions t25j,[2Z]-[29]- In particular, 
for V = —1 and cyclic b.c, we found the following values of qd—l) (i) 2 for Ly = 2 
[25], (ii) 2 for Ly = 3 [Ij, (iii) ~ 2.1548 for = 4 [2S], and (iv) ^ 2.2641 for Ly = 5 
[26]. For Ly oo, eq. formally yields gc(-l) = (3 + ^5)72 ~ 2.6180. This is 

formal since the Potts antiferromagnet is, in general, only defined for q G Z4., owing 
to the fact that the formula (II. ip can yield a negative and hence unphysical value 
for Z{G, q, v) if v is negative and q is not a non-negative integer. There are at least 
two different ways in which, as 00, the loci B could approach the limiting form 

containing qd—l): (i) the endpoints of the complex-conjugate prongs farthest to the 
right could move down and pinch the real axis and/or (ii) the point at which the locus 
B for finite Ly crosses the real axis could increase to this limiting value. These are not 
necessarily mutually exclusive; a combination of (i) and (ii) could presumably occur, 
so that the final locus would have a structure similar to that of the Ly = 2 cyclic case 
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shown in Fig. 1(a) of [25] or the Ly = 3 case shown in Fig. 7 of [26], in which several 
curves on B meet at qd—l). Alternatively, the rightmost complex-conjugate prongs 
might bend back and intersect the rest of the boundary B away from the real axis, 
forming "bubble" regions, as we found for the cyclic Ly = 4,5 strips (Figs. 8 and 9 
in [26]). 

As regards the general distribution of zeros, for a given Ly, as f increases from 
— 1 to (i.e., K = (3J increases from — cxo to 0, these zeros contract to a point 
at g = 0. This is an elementary consequence of the fact that i^' — > 0, the spin- 
spin interaction term in the Potts model Hamiltonian, Ti, vanishes, so that the sum 
over states just counts all q possible spin states independently at each vertex, and 
Z{G,q,v) approaches the value Z{G,q,Q) = g". More generally, one can inquire 
about the maximum modulus of a zero of Z{G, q, v) in the q plane as a function of v. 
It has been proved [20] that (for a graph G without loops), for the antiferromagnetic 
Potts model partition function with |1 + f | < 1, the zeros of Z{G, q, v) lie in the disk 
|g| < Cr|f I, where C ~ 7.964 and r is the maximal degree of a vertex, with r = 3 
for our honeycomb-lattice strips. Our zeros have moduli that lie considerably below 
this upper bound. For example, for v = —1, the zeros that we have calculated and 
displayed in Figs. miHlhave moduli bounded above by about 1.2 and 1.3, respectively, 
while the above-mentioned inequality would give an upper bound of |g| < 3C ~ 23.9. 

One can also plot Bq for the ferromagnetic region < v < oo. Although we have 
not included these plots here, we note that an elementary Peierls argument shows that 
the Potts ferromagnet on infinite-length, finite-width strips has no finite-temperature 
phase transition and associated magnetic long range order. Hence, for this model Bq 
does not cross the positive real q axis for < f < oo. 



6 Partition Function Zeros in the v Plane 
6.1 General 

In this section we shall present results for zeros in the f-plane, for various values of 
q, for the partition function of the Potts model on strips of the honeycomb lattice of 
widths Ly < 5 with free boundary conditions and of widths Ly = 4,6 with cylindrical 
boundary conditions. We recall the possible noncommutativity in the definition of 
the free energy for certain special integer values of q, denoted qsp, (see eqs. (2.10), 

(2.11) oim-- 

lim Urn Z{G,q,vy/'' lim lim ^(6", g, w)^/" . (6.1) 

As discussed in [3], because of this noncommutativity, the formal definition f 1 1.2 1) 
is, in general, insufficient to define the free energy / at these special points qsp, it is 
necessary to specify the order of the limits that one uses in the above equation. We de- 
note the two definitions using different orders of limits as fqn and fnq- fnq{{G}, q, v) = 
lim„^oo limq^g.p n"^ In Z{G, q, v) and fqn{{G}, q, v) = limg^g^^ lim„_oo n'^ In Z{G, q, v) 
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Figure 4: Partition-function zeros in the q plane for the Potts antiferromagnet with 
(a) V = —1.0, (b) V = —0.75, (c) v = —0.5, and (d) v = —0.25 on strips with free 
boundary conditions and several widths Ly-. 2 (□, black), 3 (o, red), 4 (+, green), 
and 5 (<0, blue), where the colors refer to the online paper. 
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(c) (d) 

Figure 5: Partition- function zeros and asymptotic loci B in the q plane for the Potts 
antiferromagnet on the he strip with width Ly — 3 and free boundary conditions, for 
(a) V = -1.0, (b) V = -0.75, (c) v = -0.5, (d) v = -0.25. 
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Figure 6: Partition-function zeros for (a) v = —1.0, (b) v = —0.75, (c) v = —0.5, and 
(d) V = —0.25 on strips with cylindrical boundary conditions and several widths Ly-. 
4 (□, black), 6 (o, red), where the colors refer to the online paper. 
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As a consequence of this noncommutativity, it follows that for the special set 
of points q = Qsp one must distinguish between (i) {By{{G},qsp))nq, the continuous 
accumulation set of the zeros of Z{G, q, v) obtained by first setting q = qsp and then 
taking n — > oo, and (ii) {Bv{{G}, qsp))qn, the continuous accumulation set of the zeros 
of Z{G,q,v) obtained by first taking n —>■ oo, and then taking q qsp. For these 
special points (cf. eq. (2.12) of |4j), 



Here this noncommutativity will be relevant for g = and q = 1- 

In Figure [7] we show the partition-functions zeros in the f-plane, for fixed values 
of q, for strips with 2 < Ly < 5 and free boundary conditions. We have displayed 
each value of g on a different plot: (a) g = 0, (b) g ~ 1, (c), g = 2, and (d) g = 3. 
The corresponding partition-function zeros for honeycomb-lattice strips with Ly=A,6 
and cylindrical boundary conditions are shown in Figure [3 Complex-temperature 
phase diagrams and associated partition function zeros were given in [1] for Ly = 2 
for free longitudinal boundary conditions. Our present calculations extend that work 
to greater strip widths. 

For our discussion we will need some results concerning the behavior of the Potts 
model on the infinite honeycomb lattice (defined via the 2D thermodynamic limit). 
The criticality condition for the g-state Potts model on this lattice is [3T]-[33] 



Since eq. (16.31) is cubic in v, it is cumbersome to write the general solution for v as 
a function of g. The following information will suffice: the equation has (i) one real 
root in v for real g < and g > 4 (ii) three degenerate real roots, f = 0, at g = 0, (iii) 
three real roots, two of which are degenerate, at g = 4: v = —2, —2, 4; and (iv) three 
distinct real roots for < g < 4. A plot of these roots is given, e.g., as Fig. 4 of [53]. In 
the interval < g < 4, the maximal root, Vhc3, which is the transition point between 
the paramagnetic (PM) and ferromagnetic (FM) phases, increases monotonically from 
Vhc3 = at g = to Vhcs = 4 at g = 4, while the middle one, Vhc2, which is 
the transition point between the paramagnetic and antiferromagnetic (AFM) phases, 
decrease monotonically from Vhc2 = at g = 0, through Vhc2 = — 1 at g = (3 + a/5)/2, 
and then to Vhc2 = — 2 at g = 4. Only the interval > Vhc2 > — 1 corresponds to 
a physical PM-AFM transition; for g = (3 + \/3)/2, this transition occurs at zero 
temperature, i.e., v = —1, and for larger values of g, the Potts antiferromagnet has 
no physical PM-AFM transition and is disordered and noncritical even at T = (e.g., 
[36]). A rigorous result is that the g-state Potts antiferromagnet on the honeycomb 
lattice is disordered and noncritical even at T = if g > 4 [37] . As noted above, 
the critical point for the Potts antiferromagnet for g ^ Z_|_ has a formal, rather than 
directly physical, significance. The lowest root, Vhci, is unphysical; it is equal to at 
g = 0, decreases through v = —2 at g = 2 to a minimum of —9/4 at g = 27/8 = 3.375 
and then increases slightly to reach the value —2 again at g = 4. In the interval 
< g < 4 it is convenient to write 



{B,{{G},qsp))nq^{B,{{G},qsp)) 



qn 



(6.2) 



— 3qv — g^ = . 



(6.3) 




(6.4) 
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with < < 7r. The main cases of interest here are contained within the discrete set 
of values q = = q{Or), where e*^'' is a certain root of unity given by 



27r 

^ - _ , r G Z , 



(6.5) 



These special values qr in were discussed by Tutte and Beraha in connection with zeros 
of chromatic polynomials [38l[39j and are also of interest since they correspond to roots 
of unity for the deformation parameter in the Temperley-Lieb algebra relevant for the 
Potts model HH SUj. Some values are = 4, 0, 1, 2, (3 + V^)/2, 3 for 1 < r < 6, 
respectively. For q = qr, the solutions of eq. (16. 3p for v have simple expressions in 
terms of trigonometric functions, which will be of use for our discussion below 



Vhci{r) = -4cos(^^j cos 



Vhc2[r 



-4 cos — 
V r 



cos 



3 Vr 

TT /I 

3f- + l 



and 



4 cos — 
\ r 



(6.6) 



(6.7) 



(6.8) 



This set is dual, via the map v ^ q/v, to the set of solutions of the corresponding 
criticality condition for the triangular lattice, + 3w^ — g = 0, viz., vt^-q = — 1 + 
2cos[2(l + ?7r)7r/(3r)] for t] = 1,0,-1, given in eqs. (27)-(29) of |12]. In previous 
studies such as [1] , it has been found that although infinite-length, finite- width strips 
are quasi-one-dimensional systems, and hence the Potts model has no physical finite- 
temperature transition for such systems, some aspects of the complex-temperature 
phase diagram have close connections with those on the (infinite) honeycomb lattice. 
We shall discuss some of these connections below. One can also express the solutions 
of eq. (16. 3p as values of q for a given v; these are 



-3v± v/^^2(4^7T9) 



(6.9) 



and are real for v > —9/4. 



6.2 q = 

From the cluster representation of Z{G, q, v), eq. (11.11) . it follows that this partition 
function has an overall factor of where k(G) denotes the number of components 
of G, i.e., an overall factor of q for a connected graph. Hence, Z{G,q = 0,v) = 0. 
In the transfer matrix formalism, this is evident from the overall factor of q coming 
from the vector w. However, if we first take the limit n ^ oo to define B for q ^ 
and then let g — > or, equivalently, extract the factor q from the left vector w, 
we obtain a nontrivial locus, namely (;B„({(j}, 0)q„. This is a consequence of the 
noncommutativity (16. 2p for q = 0. 

With the second order of limits or the equivalent removal of the factor of q in 
Z, we obtain the zeros for g = shown in Figures [7(a) (free boundary conditions) 
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Figure 7: Partition-function zeros, in the v plane, of (a) Z{G,q,v)/q for q = 0, and 
Z{G, q, v) for (b) q = 0.999, (c) q = 2, and (d) g = 3 on strips with free boundary 
conditions and several widths Ly-. 2 (□, black), 3 (o, red), 4 (+, green), and 5 (0, 
blue), where the colors refer to the online paper. 
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Figure 8: Partition-function zeros, in the v plane, of (a) Z{G,q,v)/q for q = 0, and 
Z{G,q,v) for (b) g = 0.999, (c) q = 2, and (d) g = 3 on strips with cyhndrical 
boundary conditions and several widths Ly-. 4 (□, black), 6 (o, red), where the colors 
refer to the online paper. 
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Figure 9: Partition-function zeros, in the v plane for g = gs = (l/2)(3 + y/b ) for (a) 
free boundary conditions with Ly = 2 (□, black), 3 (o, red), 4 (-f-, green), and 5 (0, 
blue); (b) cylindrical boundary conditions and Ly =: 4 (□, black), 6 (o, red), where 
the colors refer to the online paper. 
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and[8](a) (cylindrical boundary conditions). The zeros appear to converge to a roughly 
circular curve. We see in these figures that the limiting curves cross the real f-axis at 
V ~ —6. As Ly increases, the arc endpoints on the upper and lower right move toward 
the real axis. It is possible that these could pinch this axis at f = as Lj^ oo, 
corresponding to the root of (16.31) for g = 0. 

6.3 q = l 

For g = 1, the spin-spin interaction in Ti always has the Kronecker delta function 
equal to unity, and hence the Potts model partition function is given by 

^(^,5 = l,t;) = e^l^l = (l + t;)l^l (6.10) 

where \E\ is the number of edges in the graph G. This has a single zero at v = — 1. 
But again, one encounters the noncommutativity (16. 2p for q = 1. It is interesting to 
analyze this in terms of the transfer matrix formalism. At this value of q, both the 
transfer matrix and the left vector w are non-trivial. There is thus a cancellation 
of terms that yields the result (16.101) . The strip with Ly = 2 and free boundary 
conditions is the simplest one to analyze: the eigenvalues and coefficients for q = 1 
are given by 

Xi{l,v) = V* ] ci(l,0) = (6.11a) 

\2{l,v) = [l + vf- C2(l,0) = (1 + t;) (6.11b) 

Thus, only the second eigenvalue contributes to the partition function, and it gives 
the expected result Z{2 x m, FF, q = l,v) = {1 + w)^™"*"^. The strip with Ly = 3 and 
free boundary conditions is similar: there is a single eigenvalue Ai(l,f) = (1 + v)^ 
with a non-zero coefficient ci(l,f) equal to (1 + v) for odd L^ and (1 + v)^ for even 
Lrc- The other four eigenvalues including v^{l + v)^ and the roots of 

(6.12) 

have identically zero coefficients. Thus, the partition function takes the form Z{3 x 
m, FF, q = l,v) = {1 + 1))8™+'= = (i -)_ -y)!-^!^ where k=l for odd L^ and 5 for even L^. 

In general, we conclude that at g = 1, only the eigenvalue A = (1 + v^^^ con- 
tributes to the partition function for cylindrical boundary conditions, and only the 
eigenvalue A = (1 + v)^^y~^ contributes to the partition function for free b.c. The 
other eigenvalues do not contribute because they have zero coefficients in eq. (12. 2p . 
This is analogous to what we found in earlier work for cyclic strips [H [5l [11] . 

In our present case, in order to get insight into Bgn, we have computed the par- 
tition function zeros for a value of q close to 1, namely, q = 0.999 (see Figures [7](b) 
andlHl^b)). The patterns of zeros show less scatter for cylindrical, as contrasted with 
free, boundary conditions. Subsituting r = 3, i.e., g = gs = 1 in the solutions 
(I6.6p - (I6.8I) of the criticality equation (16. 3p . we have 

vj^^^ = -4 cos(7r/3) cos(27r/9) = -1.5320888... (6.13) 
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Vhc2 = -4cos(7r/3)cos(47r/9) = -0.34729635... 
Vhc3 = 4cos(7r/3) cos(7r/9) = 1.8793852... 




Our results are consistent with the inference that as L 



'y 



the real v axis at the values of Vhci and Vhc2 in eqs. There is evidently no indication 
of any crossing near the value Vhcs- A noteworthy feature of the patterns of zeros, at 
least for the case of cylindrical boundary conditions, is that they do not appear to 
exhibit the prongs that tend to occur for some other cases discussed here. 



The zeros for the q = 2 Ising case are displayed in Figures [7|(c) and [8](c) for free 
and cylindrical boundary conditions, respectively. (Fig. ^c) extends our previous 
calculation presented in Fig. 3 of [T] to greater widths.) From our earlier work [TT] 
one knows that the loci are different for strips with free or periodic transverse 
boundary conditions and free longitudinal boundary conditions, on the one hand, and 
free or periodic transverse boundary conditions and periodic (or twisted periodic) 
longitudinal boundary conditions. One anticipates, however, that in the limit of 
infinite width, the subset of the complex-temperature phase diagram that is relevant 
to real physical thermodynamics will be independent of the boundary conditions used 
to obtain the 2D thermodynamic limit. 

In the 2D thermodynamic limit, one knows the complex-temperature phase di- 
agram exactly for the q = 2 (Ising) case. (This isomorphism involves the redefini- 
tion of the spin-spin exchange constant Jpotts = 2Jising and hence -R'potts = 2Kising, 
where -ft'potts is denoted simply K here.) Since the honeycomb lattice is bipartite, the 
phase boundary separating the paramagnetic and ferromagnetic phases maps into 
that separating the paramagnetic and antiferromagnetic phases under the transfor- 
mation K —K, i.e., a —>■ 1/a, where a = = v + 1. The total boundary locus 
B is invariant under this inversion a 1/a. Because of this symmetry, it is con- 
venient to discuss the phase diagram first in terms of the variable a; the features 
in the v plane then follow in an obvious manner. Following the calculation of the 
zero- field free energy / of the Ising model on the square lattice [43j, / was calculated 
on the triangular and honeycomb lattices in [Hj. The critical points separating the 
PM and FM phases are given by Oc = 2 + y/S and l/oc = 2 — y/S, respectively. In 
terms of v, these correspond to the values Vhcs and Vhc2 with r = 4 in eqs. fl6.8p and 
fl6.7p . The complex-temperature phase diagram, with boundaries comprised by the 
locus B, was given in the plane of the variable tanh(i^r) = (a — l)/(a + 1) in (Fig. 
3 of) [15] and in the variable a in (Fig. 1(c) of) [16]. The locus B separates the 
complex a plane into three phases: (i) the physical PM phase occupying the interval 
2 — < a < 2 + v^, and its complex-temperature extension (GTE), where the 
Sq symmetry is realized explicitly {Sq being the symmetric group on q numbers, the 
symmetry group of the Hamiltonian) , (ii) the physical ferromagnetic phase occupying 
the interval 2 + < a < oo (and its CTE), and (iii) the physical antiferromagnetic 
phase occupying the interval < a < 2 — y/S and its CTE. The boundary separat- 
ing the CTE of the FM and AFM phases is an arc of the unit circle a = e*^ with 



6.4 q = 2 
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7r/3 < 9 < 57r/3; it thus has endpoints at a = e^*'^/^ and crosses the real a axis at 
a = —1. The rest of S is a closed curve crossing the real axis at a = 2±a/3 and having 
intersection points with the above-mentioned circular arc at the points at a = ±i. 
In pp the locus B for an infinite-length free and cyclic strip with width Ly = 2 were 
compared with this 2D phase diagram. 

Using our exact results, we can compare the loci for various strip widths 
and either free or periodic transverse boundary conditions with the known complex- 
temperature phase diagram for the Ising model on the infinite 2D honeycomb lattice. 
This comparison is simplest for the case of free boundary conditions, so we concen- 
trate on these results. For the finite values of Ly that we have considered, the loci 
B^ inferred from these zeros clearly contain a circular arc crossing the real v axis 
at w = —2 and have intersection points at w = — 1 ± i, just as is the case with the 
exactly known locus B for the infinite 2D honeycomb lattice. Moreover, one sees 
that as Ly increases, the endpoints of the complex-conjugate arcs move down toward 
the real axis. As Ly — > cxd, we expect that these arc endpoints will cross the real v 
axis at the points f = 1 ± that constitute the intersections, with the real v axis, 
of the complex-temperature phase boundaries B for the Ising model on the infinite 
honeycomb lattice. As in our earlier studies, this comparison shows that, although 
the behavior of the asymptotic locus B for infinite-length lattice strips is qualitatively 
different from the locus for the thermodynamic limit of the two-dimensional lattice 
as regards the physical phase transitions (owing to its quasi-one- dimensional nature), 
its features for complex temperatures show many similarities with the exactly known 
features for the 2D thermodynamic limit. 



6.5 g=(3 + v^)/2 



One of the useful features of exact solutions for Z{G, q, v) for arbitrary q and v 
is that they allow one to analyze values of q that are not positive integers and hence 
cannot be represented in Hamiltonian form, but instead via the relation fll.ip . Within 
the sequence of the q^s, the first such example is provided by the case r = 5, i.e., 
q = q5 = {3 + \/5)/2 ~ 2.6180. For this value, the criticality condition for the model 
on the (thermodynamic limit of the) honeycomb lattice has the solutions 



Vhc3 



1 + V3(5 + 2V5) 



Vhc2 



Vhcl = 77 



1 



3(5 + 2V5) 



3.165 



-2.165 



(6.16) 

(6.17) 
(6.18) 



The first of these is the PM-FM critical point, while the second formally corresponds 
to the critical temperature for the Potts antiferromagnet going to zero. The third is a 
complex-temperature singular point. In Figs. [9] we plot zeros of the partition function 
in the v plane for q = q^ and strips with free and cylindrical boundary conditions. 
From these zeros, one can infer that as Ly —>■ 00, the rightmost complex-conjugate 
arc endpoints would move in and pinch the real axis at the PM-FM value Vhcs — 3.165 
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given in eq. (16.161) . The results are also consistent with crossings on B at the other 
two points Vhc2 and Vhci as well as the values v ~ —0.6 and v ~ —4.4 which are not 
roots of the criticality condition (16.31) . 

6.6 q = 3 

In contrast to the q = 2 case, the free energy of the g-state Potts model has not 
been calculated exactly for g > 3 on any 2D (or higher-dimensional) lattice, and hence 
the corresponding complex-temperature phase diagrams are not known exactly. For 
g = 3, i.e., r = 6, the solutions of the criticality equation (16. 3p given by eqs. (I6.6p - (I6.8I) 
are 

y^^3 = 2V3cos(7r/18) = 3.4114741... (6.19) 

corresponding to the physical PM-FM phase transition point, and two other roots at 
the complex-temperature values 

Vhc2 = -2V3cos(77r/18) = -1.1847925... (6.20) 

and 

Vhci = -2v/3cos(57r/18) = -2.2266815... (6.21) 

Some discussions of the complex-temperature solutions of eq. (16. 3p and their connec- 
tions with the complex-temperature phase diagram have been given in |47j-|49j. 

The partition-function zeros in the v plane for q = 3 are displayed in Figures [3(d) 
and [8](d) for free and cylindrical boundary conditions, respectively. We expect that 
the pair of complex-conjugate endpoints in this regime will eventually converge to 
the ferromagnetic critical point Vhcs d.s Ly oo. However, obviously, an infinite- 
length strip of finite width Ly is a quasi-one- dimensional system, so the Potts model 
has no physical finite-temperature phase transition on such a strip for any finite 
Ly. The q = 3 Potts antiferromagnet is disordered on the honeycomb lattice for 
all temperatures T including T = 0, so there is no finite-temperature PM-AFM 
transition. 

In the complex-temperature interval v < —1, there are considerable finite-size 
and boundary condition effects. Because of this, in previous work, a combination of 
partition-function zeros and analyses of low-temperature series expansions was used 
[i9] ; these enable one at least to locate some points on the complex-temperature phase 
boundary. As regards the infinite 2D honeycomb lattice, because of a duality relation, 
the complete physical temperature interval < T < oo, i.e., < a < 1 of the g-state 
Potts antiferromagnet on the triangular lattice is mapped to the complex-temperature 
interval —oo < v < — g) on the honeycomb lattice (and vice versa) ^8]. Ref. [48] 
found that there is a complex-temperature singularity for the g = 3 Potts model at 
Vtri,PM-AFM,q=3 = —0.79691 ± 0.00003. From duality, the corresponding singularity 
on the honeycomb lattice is Vhc,q=z — '^/'^tri,PM—AFM,q=3 — — 3. 76454 ± 0.00015. One 
anticipates that as L^^ — > cxd for the infinite-length, width-L^ strips of the honeycomb 
lattice, the left-most arcs on will cross the real v axis at this point. There are 
several other cases of interest, such as g = 4 and g^ values with r > 7. For brevity, 
we do not consider these here. 
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7 Internal Energy and Specific Heat 



It is of interest to display some of the physical thermodynamic functions for the 
Potts model on the infinite-length limits of these strips. Having calculated the parti- 
tion function, one obtains the free energy per site, f{G,q,v) as 

f{G,q,v) = -\ogZ{G,q,v) (7.1) 
n 

for finite n, with the n — » oo limit having been defined in eq. (11. 2p above. The internal 
energy per site, E, is 

^(G,g,^^) = -|| = -J(^^ + i)f (7.2) 



and the specific heat per site, C, is 

dE 



df ^ ^ d'f 
dv dv"^ 



(7.3) 



As the strip width Ly — ^ oo, these approach the internal energy and specific heat for 
the infinite 2D honeycomb lattice. For convenience we define a dimensionless internal 
energy 

^' = -7 = '" + ^)|- f^-^' 

Note that sgn(£'r) is (i) opposite to sgn(ii^) in the ferromagnetic case where J > 
for which the physical region is < f < oo and (ii) the same as sgn(ii^) in the 
antiferromagnet case J < for which the physical region is — 1 < u < 0. Of 
course, the infinite-length limits of the honeycomb-strips considered here are quasi- 
one-dimensional systems, so that /, E, and C are analytic functions of temperature 
for all finite temperatures. 

We recall the high-temperature (equivalently, small-|if|) expansion for an infinite 
lattice of dimensionality d > 2 with coordination number A: 



J 2 



q q'^ 



(7.5) 



Here, A = 3 for the infinite honeycomb lattice. Again, we recall that in papers on the 
q = 2 Ising special case, the Hamiltonian is usually defined as Hi = —JiJ2{ij) ^i^j 
with (Tj = ±1 rather than the Potts model definition of H, so that one has the 
rescaling 2Ki = K, where Ki = (3Ji. Furthermore, Ei = —J{aiaj) raher than 
the Potts definition E = —J{Saiaj), where (ij) are adjacent vertices. Hence, for 
example, for q = 2, with the usual Ising model definitions, Ei{v = 0) = rather than 
E = — JA/(2g) and the high-temperature expansion is Ei = —J{A/2)[K + 0{K^)] 
rather than the q = 2 form of fl7.5p . Similarly, for T — > 0, with the conventional 
Ising definition, Ei —\J\ for both the ferromagnetic and antiferromagnetic cases, 
while with our Potts-based definition, E —(A/2) J = —3 J/2, i.e., Er 3/2 for 
the ferromagnet and E ^ for the antiferromagnet. 
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We next present plots of the the (reduced) internal energy Er per site and the 
specific heat per site C on infinite-length honeycomb-lattice strips for three values 
of q in increasing order, q = 2, q = {3 + v^)/2, and g = 3 in the respective Figs. 
[TUlfT^ Each plot contains curves for Ly from 2 to 5 for strips with free boundary 
conditions and for Ly = 4 and Ly = 6 for strips with cylindrical boundary conditions. 
The free energy and its derivatives with respect to the temperature are independent 
of the longitudinal boundary conditions in the limit L^, — > cxd, although they depend 
on the transverse boundary conditions |1] . As expected, in the vicinity of the infinite- 
temperature point V = 0, the results for the internal energy are well described by the 
first few terms of the high-temperature series expansion given in eq. (17. 5p . One can 
see the approach of Er to its zero-temperature limit of 1.5 for the ferro magnet as v 
increases through positive values. This approach is more rapid for the strips with 
cylindrical boundary conditions, as is understandable since these minimize finite-size 
effects in the transverse direction. One can also see the approach of Er to its zero- 
temperature limit of zero for the antiferromagnet as v decreases toward v = —1. 

With regard to the specific heat, the plots show maxima which occur at values of 
V, denoted Vm, that depend on the values of q and Ly and the transverse boundary 
conditions. In the ferromagnetic case, these approach the critical values for the infinite 
honeycomb lattice, vpm-fm, as Ly increases. For the results shown, this approach is 
from above (below) in v for the case of free (cylindrical) boundary conditions. The 
heights of the maxima increase as Ly increases, in accordance with the fact that on 
2D lattices, the specific heat diverges at the PM-FM critical point for the values 
of q shown. (This divergence is logarithmic for q = 2 [43]; more generally, for the 
interval < g < 4 where the 2D Potts ferromagnet has a second-order transition, 
the specific heat exponent is given by a = a' = (2/3) (tt — 29) /{n — 6) [231 El], where 
6 = 2 arccos(g^/^/2) as in eq. (16.41) . so a = 2/9 for g = (3 + V5)/2 and a = 1/3 for 
g = 3.) This behavior as a function of increasing strip width Ly is the analogue, for 
these infinite-length strips, of the standard finite-size scaling behavior of the specific 
heat on L X L sections of a regular lattice, for which \vm — vpm-fm\ ~ L~^^'^ and 
C{v = Vrn) ~ L"^'^, where u = u' is the correlation length critical exponent [501 \5T] 
for the Potts ferromagnet on a two-dimensional lattice, with being related to a by 
the hyperscaling relation du = 2 — a, so that u = 1 for g = 2 and u = 5/6 for g = 3. 
Note that a < for < g < 2, so that, although the Potts ferromagnet has a PM-FM 
critical point for this range of g, with a divergent correlation length, the specific heat 
has only a finite, rather than divergent, nonanalyticity at the critical point. 

We next consider the plots of the specific heat in the case of the Potts antiferro- 
magnet. For q = 2, C exhibits a maximum at a value of v that approaches the value 
vpM-AFM = 1 — \/3 for the 2D lattice as Ly increases, and the height of the maximum 
increases; these are simply related, by — > —K, to the behavior for the g = 2 ferro- 
mag net. For q= {3 + ^/l)/2 the curves for the specific heat in the antiferromagnetic 
region of Fig. [11] exhibit a maximum not too far from the zero-temperature value 
V = —1. However, the height of the maximum does not increase very much as Ly 
increases, for either the free or cylindrical boundary conditions. Finally, we show Er 
and C for g = 3, in Fig. [12] For this value of g the Potts antiferromagnet has no 
finite-temperature PM-AFM transition and is noncritical even at T = 0; consistent 



30 



with this, although C has a maximum, the height of this maximum does not increase 
with increasing Ly. 




Figure 10: Reduced internal energy = —E/J and specific heat as functions 

of the temperature-like variable v for the q = 2 Potts model on the honeycomb- 
lattice strips of width 2 < Ly < 5 with free boundary conditions (a) (c) and of width 
Ly = 4,6 with cylindrical boundary conditions (b) (d). The four curves shown for the 
case of free boundary conditions correspond to Ly = 2,3,4,5 as one moves upward 
and the two curves for cylindrical boundary conditions correspond to Ly = 4,6 as 
one moves upward. The plot includes both the ferromagnetic and antiferromagnetic 
Potts models, for which the temperature ranges are < f < oo and — 1 < < 0, 
respectively. 
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Figure 11: Reduced internal energy = —E/J and specific heat C/ks as functions 
of the temperature-hke variable v for the Potts model with q = {S + V5)/2 on the 
honeycomb-lattice strips of width 2 < Ly < 5 with free boundary conditions (a) (c) 
and of width Ly = 4,6 with cylindrical boundary conditions (b) (d) . Ordering of 
curves is as in Fig. [TOl 
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(a) 



(b) 




Figure 12: Reduced internal energy Er = —E/J and specific heat C//cb as functions 
of the temper at ure-hke variable v for the g = 3 Potts model on the honeycomb- 
lattice strips of width 2 < Ly < 5 with free boundary conditions (a) (c) and of width 
Ly = 4, 6 with cylindrical boundary conditions (b) (d). Ordering of curves is as in 

Fig. m 
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8 Appendix: Transfer Matrix for Ly = A Strip with 
Cylindrical Boundary Conditions 

The number of elements in the basis is eight: P = 2,3,4, ^1,2,3 + ^1,2,4 + (^1,3,4 + 

(52,3,4, 51,2^3,4, ^1,2 + ^3,4, ^1,3 + '^2,4, ^1,4^2,3, ^2,3 + ^1,4, !}• The transfer matrix is given 
by 

T = 



/ 




iv''DiD3S25 




2v^DiTi4 


2v^Sl, 


v^D2Dl 


2v^D3S23 




\ 




v''DlD3S25 


V^DiT22 


v^DlT2:i 


v^DiT24 


V^S25T25 


V^D2D3S25 


v'^T27 


f^7'28 














2v^Di 



















2v^DiS25T25 




v'^DiTii 




v'D2Si, 


21^5^25747 










2v^DiS25T25 




2v^D^T^i 




v'D2Si, 


21^5^25747 










^v'^DiTq2 


v'^DlFi 


2v^DiF2Si2 


2v^F2Te5 




2v''Te7 


v'Sl2 






v^DfTn 


2v^DiTt2 




2v^DiT74, 


2v^Tr5 




v^Tjr 






v 


DIT23T81 


4DiTs2 




2DiTsi 


2^85785 


Ts6 


787 


788 


/ 



(8.1) 



where the factors D^. and are defined in eqs. fl3.8ap and (13. Sb^ : the Tij are given 
by 

Tig = 2q + 6v + v^ (8.2a) 

Ti4 = + 4qv + 7v^ + (8.2b) 

Ti5 = q + 8v + Av^ (8.2c) 

Ti8 = 2q'^ + 6qv + 8v'^ + (8.2d) 



^22 


= 5q^ + q^v + Mqv + lOgw^ + 65w^ + 33w^ + Av^ 


(8.3a) 


^23 


= q^ + 5qv + 8f ^ + 


(8.3b) 


^24 


= q^ + 7q\ + 17qv'^ + 19v'^ + 2v^ 


(8.3c) 


T25 


= q'^ + 6qv + llv"^ + 2v^ 


(8.3d) 


^27 


= 5q^ + q\ + 28qv + 8qv'^ + AW + 19w^ + 2v^ 


(8.3e) 


^28 


= q^ + 6q\ + Uqv'^ + llv'^ + V* 


(8.3f) 



T43 = q^ + Gq'^v + Uqv^ + Uv^ + (8.4a) 

T44 = q^ + 8q\ + 26q^v^ + A2qv^ + 33v^ + 2v^ (8.4b) 

T47 = q^ + 5qv + 7v^ + (8.4c) 

T48 = q'^ + 7q^v + 20q^v^ + 28qv^ + 19v^ + v^ (8.4d) 



T53 = 2q'^ + 8qv + lOv^ + (8.5a) 

T54 = q^ + Gq'^v + 13qv^ + 12v^ + v'^ (8.5b) 

T55 = + 8qh + 32qh^ + 2q%^ + 68qv^ + 12gw^ + Qlv"^ + 22v^ + 2v^ (8.5c) 

T58 = 2g3 + 10g\ + 18gf 2 + 14f 3 + t;^ (8.5d) 
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^61 


= q + 6v + 2v'^ 


(8.6a) 


262 


^ + 6qv + qv'^ + lOv'^ + 3v^ 


(8.6b) 


265 


= 5^ + Aqv + 8v'^ + 2v^ 


(8.6c) 


Tm 


^ q^ + Sqv + 3qv^ + 2lv^ + 16^;^ + ?>v^ 


(8.6d) 




= + Iq^v + q^v'^ + I9>qv'^ + Aqv^ + llv^ + 


(8.6e) 



Tyi = + 8g^^; + gV + 25?^;^ + 6?^;^ + 32^;^ + 13^;^ + t;^ (8.7a) 

= 25^ + 17?^T; + gV + 62gV + 8gV + 117?T;^ + 25gi;^ + 98T;^ + 34v^ + 2t;6 

(8.7b) 

T73 = g^ + 7g^t; + 21gV + 32gf''^ + 22f^ + t;^ (8.7c) 

T74 = + Sg^t; + 28gV + 5Vf ^ + 58gt;^ + 30?;^ + (8.7d) 

T75 = g^ + 9gS + 37gV + gV + 87gV + 7gV + 118A + 19^;^? + 74^;^ + 21f^ 
W (8.7e) 

Tye = q* + IQvq^ + 2^^^^ + 4:lv'^q^ + 16w^g^ + w^g^ + 88w^g + 52t; ^ + 7gt;-' + 88t;^ 
+72^;^ + llv^ + v'^ (8.7f) 

7^7 = + 19gS + q^v^ + 78g^i;^ + 8g^i;^ + 180?^^;^ + 2%q\^ + 240?^;^ + 52^1;^ 

+149^;^ + 47^;^ + 2v'^ (8.7g) 
T78 = g^ + 9v?^ + 36?V + 83gV + 118?V + 99?^;^ + 41^;^ + 't;^ (8.7h) 



Tgi = q^ + bq% + Uqv"^ + qv^ + IQv^ + Av* (8.8a) 

= g^ + llg^f + 55gV + gV + 162gV + 9gV + 302gV + 34gV + 345gw^ 
+67gt;^ + gf^+ 192f^ + 61f^ + 4f^ (8.8b) 

Tga = + lOg^t; + Abq^v'^ + llSg^t;^ + lMq%^ + 194gt;^ + gt;^ + 96t;^ + Av'^ (8.8c) 

= + llg% + SSg^-u^ + 163g^'t;^ + 313g^'t;^ + 396g2't;^ + 313gv'^ + qv'^ + 126v^ 
+Av^ (8.8d) 

Tss = q^ + 6q^v + Uqv^ + 13v^ + t;^ (8.8e) 

= g^ + 6g^w + 16gV + 25gv^ + gv^ + 22v^ + 4^^ (8.8f) 

Tge = g^ + 12g^'i; + qh^ + 65gS^ + 12gS^ + 206g^v^ + 63g2'i;^ + 2g3'i;^ + 412g^^;^ 
+184g2-t;5 + I5q%^ + 512qv^ + 312qv^ + 46g^;^ + g^;^ + 320v^ + 256v^ 
+60v^ + iv^ (8.8g) 

Tsr = 2g^ + 24g% + 132g^f ^ + 2f ^g^ + 434g\^ + 18g%^ + 934g^t;^ + 72g''^t;^ 
+ 1346g2't;^ + 166g2't;^ + 1232gv6 + 230gv^ + 2qv^ + 560v^ + 160v^ + 8v^ 

(8.8h) 

Tgg = g^ + 12g^'i; + 66q%^ + 218g^'i;^ + 477gS^ + 718g^v^ + 739g^v^ + 486gv^ 

+g^;^ + 165^;^ + Av^ (8.8i) 
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The vectors w and Uid are given by 



w, 



odd 



Uid = 



where 



q {Dl ^D^F^, qDl 2qD^F^, 2Ff , q + 2v + v\ 2Ff, qF^) (8.9a) 

q {DlXl DlX^, 2D,X,, 2Xf, {q + 2v + v^)Xl 2X2X,, X,) 

(8.9b) 

(1,0,0,0,0,0,0,0) (8.9c) 



^3 



+ 6v^ + IBv'^q + 20v^q'^ + 15v^q^ + 6gS + q^ (8.10a) 

v'^ + 7v^ + 21v^q + 35v^q'^ + 35g\^ + 21q\'^ + 7vq^ + q^ (8.10b) 

+ 5v^ + lOv^q + lOv^q'^ + 5vq^ + (8.10c) 

+ 8v^ + 28v^q + 56v^q^ + 70q\^ + 56q%^ + 28v'^q^ + 8q% + q' 

(8.10d) 
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